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^' Abstract 

We compute the low lying spectrum of the overlap Dirac operator in the 
r> I deconfined phase of finite-temperature quenched gauge theory. It suggests 

C^ ' the existence of a chiral condensate which we confirm with a direct stochastic 

estimate. We show that the part of the spectrum responsible for the chiral 
condensate can be understood as arising from a dilute gas of instantons and 
anti-instantons. 
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I. INTRODUCTION 

The fermion spectrum near zero eigenvalue is closely related to gauge field topology and to 
cliiral symmetry breaking. While there is considerable experimental and theoretical support 
for chiral symmetry breaking in gauge theories with dynamical quarks at zero temperature 
and theoretical arguments for its restoration above a critical temperature, the situation is 
less clear for the nominally simpler case of fermions in the background of quenched gauge 
fields. To improve our understanding of the quenched, deconfined phase, we have studied 
the spectrum of the hermitian overlap Dirac operator |jl| in that region. We find a segment 
of the spectrum concentrated at and near zero eigenvalue and separated from the bulk of the 
spectrum. The bulk of the spectrum begins to rise rapidly at larger eigenvalues. The exactly 
zero eigenvalues are associated with the global topology of the gauge field configurations. 
The statistical properties of the small eigenvalues are in correspondence with predictions 
from a dilute gas of instantons and anti-instantons. Small, non-zero eigenvalues with these 
properties give rise to a finite chiral condensate. 

Gauge field topology plays a central role in QCD. The presence of gauge field configura- 
tions with non-trivial topology indicates that massless fermions will have exact zero modes. 
These zero modes cause an explicit breaking of the axial U(l) symmetry and result in a 
massive i]' [^. Conventional wisdom says that the axial U(l) symmetry remains broken at 
all temperatures since one does not expect a complete suppression of non-trivial gauge field 
backgrounds. Recent studies using dynamical staggered fermions indicate that the axial 
symmetry most likely remains broken at high temperatures although the magnitude might 
be considerably smaller than at low temperatures ^^. Since topology is also expected to 
be suppressed at high temperatures 0, this result is consistent with expectations. Chiral 
symmetry breaking, on the other hand, comes from the finite density of eigenvalues near 
zero 1^, so we expect this density to be zero at high temperature in full QCD since chiral 
symmetry is restored in that case. 

The most likely scenario for the spectrum of the massless Dirac operator at high temper- 
atures is to have a delta function at zero due to topology, followed by a gap and a continuous 
spectrum of eigenvalues resulting in a theory with unbroken chiral symmetry and a broken 
axial U(l) symmetry. If we adopt the instanton picture for topology, we expect a dilute gas 
of instantons and anti-instantons at high temperature since the topological susceptibility is 
highly suppressed. However, this dilute gas of instantons and anti-instantons should not 
give rise to a chiral condensate in high-temperature, full QCD. A natural explanation in 
the context of instanton models is that instantons and anti-instantons form molecules at 
high temperatures P]. One expects that this formation of molecules is primarily due to 
interactions induced by fermions [§|,3. 

In this paper, we will study the spectrum of the massless overlap Dirac operator [Q on 
the lattice in pure SU(2) and SU(3) gauge theories on lattices with A''^ = 4 in the deconfined 
phase. We will study several different ensembles to understand the finite volume effects and 
also the effect of going deeper into the deconfined phase. In addition to directly studying 
the spectrum, we will also study the chiral condensate and scalar susceptibility. In some 
ensembles, we will also compare our results with those obtained using staggered fermions 
under these conditions. 

We have the following results: 



• The overlap Dirac operator has exact zero eigenvalues indicating that gauge field con- 
figurations with non-trivial topology persist in the high temperature phase. 

• The spectrum of the non-zero eigenvalues of the overlap Dirac operator has two parts 
separated by a region in eigenvalue where the density is essentially zero. 

— The distribution of very small non-zero eigenvalues of the overlap Dirac operator 
is well-described by a dilute gas of non-interacting instantons and anti-instantons. 

— The bulk part of the spectrum begins at larger eigenvalue and rises smoothly and 
steeply. 

• The very small eigenvalues of the overlap Dirac operator have the properties that are 
needed to give a chiral condensate. The separation between the very small eigenvalues 
and the bulk of the spectrum results in a minimum in the scalar susceptibility. 

These features are quahtatively different from the staggered Dirac operator in a quenched 
theory. The staggered Dirac operator does not have any exact zero eigenvalues. Studies with 
dynamical staggered fermions have indicated that chiral symmetry is restored at high enough 
temperatures in QCD with Nf flavors of massless quarks and also in the quenched theory. 
Early numerical studies of chiral symmetry restoration at high temperatures using staggered 



fermions in the "quenched" approximation can be found in Ref. |I0|. Further studies using 
dynamical staggered fermions have established a phase transition showing chiral symmetry 
restoration at high temperatures. Recent reviews of the lattice results can be found in 



Ref. [|TT|, and Ref. |T^ contains the most recent results for the Nt = 4 phase transition. 

The lack of exact zero eigenvalues in the staggered fermion spectrum follows from the 
breaking of the continuum chiral and flavor symmetry at flnite lattice spacing |TB|. The 



observed restoration of chiral symmetry at high temperature is consistent with the prediction 
of Ref. P] for SU(2) gauge theory with staggered fermions. However, a gap in the spectrum 



of the staggered Dirac operator at high temperature has not been convincingly established, 
see e.g. Ref. |^. A tail of small eigenvalues seems to persist, possibly consisting of the 
"shifted" would be zero modes due to global topology. 



The proof of chiral symmetry restoration at high temperature in Ref. [Q is for dynamical 
fermions, while we are only studying the quenched theory here. So, we could explain the 
presence, and even accumulation, of very small eigenvalues in the spectrum of the overlap 
Dirac operator as a quenched artifact. This is not satisfactory, however, since staggered 
fermions in a quenched theory show no such accumulation of small eigenvalues, but rather 
a tapering tail of small eigenvalues, indicative of a chirally symmetric phase. As is well 
known, massless fermions are non-trivial to construct on the lattice [0]. The overlap Dirac 
operator is one solution to the problem. It is not ultra-local, i.e. its interaction size is not 
flnite, but its interactions are exponentially decreasing with distance. It therefore does not 



fit into the assumptions used in Ref. [0] , and the proof presented there need not hold even 
in the dynamical theory with overlap fermions. We do not address the issue of the existence 
of a chirally symmetric phase in the dynamical theory with overlap Dirac fermions in this 
paper. 

Care has to be taken in a finite temperature study in the quenched approximation. In 
the pure gauge case, gauge field configurations with the Polyakov loop phase near any Z^r 



phase are equivalent, related by a global Zn symmetry. However, the coupling to fermions 
is not Z]sr symmetric and the phase with the Polyakov loop along the positive real axis is 
preferred. Configurations with "non-trivial expectation value" for Polyakov loop phase can 
give an "unphysical" signal for a chiral condensate 0] (negative for SU(2) and complex for 
SU(3) with anti-periodic boundary conditions in the time direction for fermions). In pure 
SU(N) gauge theories, configurations with expectation values for the Polyakov lines related 
by a simple Zn factor are equally likely. We will therefore generate pure gauge ensembles so 
that every configuration has a positive expectation value for the Polyakov loop and impose 
anti-periodic boundary conditions in the time direction for fermions. So in a sense, this is 
not a strictly quenched calculation, but has included in it by hand one important feature of 
the theory with dynamical quarks. 

The organization of the paper is as follows. In the second section, we will briefiy sum- 
marize what we need to know about the overlap Dirac operator. We will present our results 
for the eigenvalue spectrum in Section III and focus on the very small eigenvalues of the 
overlap Dirac operator. We will show that the spectrum of the very small non-zero and ex- 
act zero eigenvalues is described quite well by a non-interacting dilute gas of instantons and 
anti-instantons. We use the word "dilute" to emphasize that topology is highly suppressed 
and that the spectrum of these very small eigenvalues is separated from the bulk of the 
spectrum. We will present our data for the chiral condensate and scalar susceptibility and 
discuss the relationship between the small, non-zero eigenvalues and the chiral condensate 
in Section IV. Our conclusions are in Section V. 

Similar studies of the fate of chiral symmetry breaking in the deconfined phase of 



quenched QCD have been done with domain wall fermions by the Columbia group p6[ 
who measured (ipip) and found evidence of topology from its increase at small quark mass 
and by Lagae and Sinclair [^ who studied low eigenvalues and meson propagators and 
found evidence of topology from both. 



II. OVERLAP DIRAC OPERATOR 



The massive overlap Dirac operator is given by |1|,|T8 

/}(/.) = i[l + /.+ (l-/i)75e(if J] (1) 

with < ;U < 1 describing fermions with positive mass all the way from zero to infinity. 
The hermitian operator H^ is just 75-D^„ where D^ is the usual Wilson-Dirac operator with 
a negative mass on the lattice under consideration. 



The propagator for external fermions is given by [0,|l| 



D-\f,) = {1 - f,)-' [D-\f,) - 1\ . (2) 

In many cases, it is more convenient to use the hermitian version: Ho{fi) = 75D(yu). It 
is easy to see that 

H'oif^) = (i - f^')H'M + f^'; mo),i,] = 0. (3) 



Each eigenvalue < A^ < 1 of H^{0) is doubly degenerate with opposite chirality eigenvec- 
tors. In this basis, Ho{fi) and D{fi) are block diagonal with 2x2 blocks, e.g., D{fi): 



(l-/i)A2 + /i (l-^)AVr^A^^ 



-(l-^)Av^r^A2 (l-;i)A2 + ;i 

where 

For a gauge field with topological charge Q ^ 0, there are, in addition, \Q\ exact zero 
modes with chirality sign((5) paired with eigenvectors of opposite chirality and eigenvalue 
equal to unity. These are also eigenvectors of Ho{fJ.) and D{fi): 

^(/^)zcro sector : ( Q 1 ) °^ ( Q ^ j ^^^ 

depending on the sign of Q. 

In the chiral eigenbasis of iJ^(O), the external propagator takes the block diagonal form 
with 2x2 blocks 



and in topologically non-trivial background fields, the \Q\ additional blocks are 

depending on the sign of Q. 

The non-topological contribution to the fermion bilinear, which is the part that will 
survive the infinite volume limit, is 

J , 1 ^ 2^.(1 - A^) , 



A>0 



The chiral condensate is dominated by the small, non-zero eigenvalues. In the thermody- 
namic limit, it is given by the density of eigenvalues at zero p{0~^). 

A physical quantity that is more sensitive at small quark masses than the chiral conden- 
sate is the connected part of the scalar susceptibility. The non-topological contribution to 
this quantity is 

This quantity is particularly sensitive to the fi dependence of the cancellations between the 
A^ and /i^ terms in the second factor of the numerator. 



TABLE I. Pure gauge ensembles studied. All ensembles have an extent of Nt = 4 in the 
short direction. The number of configurations, the linear extent in the spatial directions and the 
estimated temperature are also given. 



Gauge group 


(3 


N 


T/T, 


L 


SU(2) 


2.3 


50 


1.0 


8 


SU(2) 


2.4 


100 


1.4 


8 


SU(2) 


2.4 


200 


1.4 


16 


SU(2) 


2.5 


200 


2.0 


16 


SU(3) 


5.71 


200 


1.03 


16 


SU(3) 


5.75 


400 


1.13 


8 


SU(3) 


5.75 


200 


1.13 


12 


SU(3) 


5.75 


400 


1.13 


16 


SU(3) 


5.85 


400 


1.38 


8 


SU(3) 


5.85 


206 


1.38 


16 



Since chiral symmetry is exact on the lattice for the overlap Dirac operator, the pion 



susceptibility is Xtt = - < i^i^ > P^) and we can define a quantity 



iU 



Xtt X 



10 






vxEimi 



f^\l - Af) 



2\2 



/i2)+/i2]2 



(11) 



This is a measure of the U{1)a symmetry breaking in a dynamical theory. It is strictly 
positive and is not sensitive to cancellations like those within the sum for Xao in Eq. (|10|). 
Due to the higher powers of the quark mass and eigenvalues in the sum, it is more sensitive 
than the chiral condensate to small eigenvalues. 

III. EIGENVALUE SPECTRUM 



We will investigate the quark spectrum on the ensembles of pure SU(2) and SU(3) gauge 
theory configurations listed in Table |. At Nt = 4, the critical coupling is (3c = 2.2986(6) 
for SU(2) and /3c = 5.6925(2) for SU(3) pf. All ensembles in Table | are in the deconfined 
phase. 

The low lying eigenvalues of Hq (including the exact zero modes) and DsD\ {Dg is the 
massless staggered operator) were computed using the Ritz variational technique [^. In 
order to deal numerically with the overlap Dirac operator, it is necessary to use a represen- 
tation for e{Hw), and we used the optimal rational approximation [^. We compare the low 
lying spectrum of the staggered and overlap Dirac operator for two typical ensembles, SU(2) 
at /3 = 2.4 and SU(3) at /5 = 5.75 both on 16^ x 4 lattices, in Fig. |l|. We see a remarkable 
difference. Both staggered and overlap spectrum show a rapidly rising bulk. But while the 
staggered spectrum has a rapidly decreasing tail of small eigenvalues, the overlap spectrum 
shows an accumulation of very small, but non-zero, eigenvalues. 

The spectral distributions of the overlap Dirac operator for all the ensembles in Table | 
have the same qualitative features: 
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FIG. 1. Histogram of the low lying eigenvalues of the staggered and overlap Dirac operator 
for the SU(2) ensemble at /? = 2.4 and SU(3) ensemble at /? = 5.75. Both lattices have a spatial 
volume of 16^. For the overlap Dirac operator the exact zero modes are not shown. 
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• There is a typical scale (A = 0.05) independent of the ensemble that separates the 
spectrum into two parts: "small" (< 0.05) and "large" (> 0.05) eigenvalues.^ 

• Small eigenvalues occur in all topological sectors. 

• The number of small eigenvalues per unit lattice volume remains roughly constant at 
fixed coupling under changes of the volume and decreases as one goes deeper into the 
deconfined phase. 

In this section, we focus on the statistical properties of the low-lying eigenvalues of the 
overlap Dirac operator. Our conclusion is that it is roughly consistent with the data to 
associate each nonzero pair of small eigenvalues with an instanton — anti-instanton pair. 
Topics related to the magnitude of the eigenvalues are discussed in the next section. 

We approach the hypothesis in steps of increasing specificity. Let n+ and n_ be the num- 
ber of separated, localized instantons and anti-instantons respectively. Then the topological 
charge of the configuration, as determined by the fermions, is Q = r2_|_ — r2_ and the total 
number of objects is n = n+ -f n_ . 

For infinite separation between instantons and anti-instantons, this gives n zero eigen- 
values. At large but finite separation, we expect \Q\ exactly zero eigenvalues and n—\Q\ = 
2 min(?T,_|_,n_) small eigenvalues. Since the nonzero eigenvalues come in ±A pairs, we as- 
sociate each nonzero pair of small eigenvalues with an instanton — anti-instanton pair. The 
number of exact zero modes of Hq in a fixed gauge field background is equal to the net 



number of levels crossing zero in the spectral flow of H^u p2[ . In addition, individual levels 
crossing zero in the spectral flow of H^ can be associated with instantons in the background 
gauge field p3| , |2^ . We performed a spectral flow on several of the ensembles in Table |I| and 



found that there is essentially a one-to-one agreement: If n^{n^) levels crossed zero from 
above (below), then there were n — \Q\ = 2 min(n+,n_) small eigenvalues for the overlap 
Dirac operator in addition to the rigorously expected |Q| = \n^ — n_ \ exact zero eigenvalues. 
Therefore our association of small eigenvalues of Hq with instanton — anti-instanton pairs is 
well justified. 

For finite temperature with T > Tc, the distribution of instanton sizes increases steeply 
with instanton size until it is cutoff on the large size end by Nt- Thus most of the instantons 
have a size not too far from Nt- We will always use a simple model in which the instantons 
and anti-instantons are assumed to be at sufficiently large separation that their interactions 
due to the pure gauge action can be neglected. 

With the assumption that the n objects are noninteracting, the number per configuration 
should have a Poisson distribution 

P{n, (n)) = (n)"e-<"V^! (12) 

where (n) is the average of the distribution. For a Poisson distribution, the average and 
variance are equal. We present our data for the distribution of n as counted by the fermions 



^ We note here that the smallest eigenvalue for the free overlap Dirac operator with anti-periodic 
boundary conditions in the temperature direction is 0.238 and 0.200 for a Wilson mass of 1.5 and 
1.8 respectively at Nt = 4. 
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TABLE II. Distribution of the number of topological objects along with predictions from a 
Poisson distribution. 





SU(2) 


SU(2) 


SU(2) 


SU(2) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


p 


2.3 


2.4 


2.4 


2.5 


5.71 


5.75 


5.75 


5.75 


5.85 


5.85 


L 


8 


8 


16 


16 


16 


8 


12 


16 


8 


16 


N 


50 


100 


200 


200 


200 


400 


200 


400 


400 


206 


n = 


11(10) 


74(75) 


27(28) 


130(130) 


1(1) 


286(290) 


83(77) 


38(35) 


379(378) 


144(137) 


n = l 


13(16) 


23(22) 


57(55) 


57(56) 


5(6) 


100(93) 


63(73) 


76(85) 


19(22) 


43(56) 


n = 2 


14(13) 


3(3) 


58(54) 


10(12) 


8(16) 


13(15) 


39(35) 


105(104) 


2(1) 


17(11) 


n = 3 


8(7) 




32(35) 


3(2) 


28(27) 


1(2) 


11(11) 


95(84) 




1(2) 


n = 4 


2(3) 




16(17) 




46(35) 




3(3) 


46(51) 




1(0) 


n = 5 


2(1) 




4(7) 




38(35) 




1(1) 


25(25) 






n = 6 






3(2) 




34(30) 






9(10) 






n = 7 






3(1) 




17(21) 






5(4) 






n = 8 










10(14) 






1(1) 






n = 9 










6(8) 












ri= 10 










2(4) 












n= 11 










3(2) 












n= 12 










2(1) 












(n)/y 


1.66 


0.29 


0.25 


0.05 


0.63 


0.32 


0.28 


0.31 


0.06 


0.05 


Ratio 


0.97 


1.09 


0.93 


0.99 


1.15 


1.08 


0.92 


1.02 


0.90 


0.83 



in Table ||. The average of the distribution per 8^ lattice volume is shown in the same table. 
The numbers in parenthesis are the predictions of the distribution using the average value 
in (|12]). We also write down the ratio of the variance and the average in the table. The data 
are in fairly good agreement with the Poisson distribution. 

To the extent that the n zero and small eigenvalues can be associated with any sort of 
objects localized in 3-space or Euclidean 4-space, the number of objects and n should, at 
fixed Nt and fixed coupling g'^, be proportional to the spatial volume V in lattice units. 
This is roughly supported by the SU(2) and SU(3) data in Table |I[ If the localized objects 
are instantons, then we also expect the density to decrease rapidly with decreasing g'^. This 
is also seen in the data. 

At the next level of detail, we can distinguish between the contributions to n from n_|_ and 
n_. If the instantons and anti-instantons are thrown into the configurations independently 
(without interaction), then for a fixed n the relative probabilities in the (n_, n+) distribution 
will be given by the binomial coefficients 



i?(n+,n_|?T,) 



n\ 



2"ni !n_! 



(13) 



Combining this with (n) values from the data gives predictions for the numbers n+ and n_ 
of positive and negative chirality states. We combine (n+,n_) with (n_,n+) and compare 
the data and predictions in Table |ni| . In almost all of the cases, the difference is within the 
statistical error expected from the value of the predicted number. 
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TABLE III. Distribution of the number of instantons for a fixed number of topological objects 
( less than eight ) along with predictions from Binomial distribution using the Poisson average. 





SU(2) 


SU(2) 


SU(2) 


SU(2) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


/3 


2.3 


2.4 


2.4 


2.5 


5.71 


5.75 


5.75 


5.75 


5.85 


5.85 


L 


8 


8 


16 


16 


16 


8 


12 


16 


8 


16 


N 


50 


200 


200 


200 


200 


400 


200 


400 


400 


206 


n = 2,n+ =0 


6(7) 


2(2) 


33(27) 


6(6) 


1(8) 


6(8) 


20(18) 


54(52) 


2(0) 


8(6) 


n = 2,n+ = 1 


8(7) 


1(2) 


25(27) 


4(6) 


7(8) 


7(8) 


19(18) 


51(52) 


0(0) 


9(6) 


n = 3, n+ = 


1(2) 




12(9) 


0(0) 


7(7) 


0(0) 


3(3) 


25(21) 




0(0) 


n = 3, n+ = 1 


7(5) 




20(27) 


3(1) 


21(21) 


1(1) 


8(8) 


70(63) 




1(1) 


n = 4,n+ = 


1(0) 




2(2) 




9(4) 




0(0) 


9(6) 




0(0) 


n = 4, n+ = 1 


1(2) 




6(9) 




18(17) 




2(1) 


30(26) 




1(0) 


n = 4, n+ = 2 


0(1) 




8(6) 




19(13) 




1(1) 


7(19) 




0(0) 


n = 5,n+ = 


0(0) 




0(0) 




1(2) 




0(0) 


2(2) 






n = 5,n+ = 1 


1(0) 




1(2) 




10(11) 




0(0) 


9(8) 






n = 5,n+ = 2 


1(1) 




3(4) 




27(22) 




1(0) 


14(16) 






n = 6, n+ = 






0(0) 




1(1) 






0(0) 






n = 6, n+ = 1 






0(0) 




7(6) 






3(2) 






n = 6, n+ = 2 






2(1) 




19(14) 






2(5) 






n = 6, n+ = 3 






1(1) 




7(9) 






4(3) 






n = 7, n+ = 






0(0) 




1(0) 






0(0) 






n = 7, n+ = 1 






0(0) 




1(2) 






0(0) 






n = 7, n+ = 2 






1(0) 




3(7) 






1(1) 






n = 7, n+ = 3 






2(0) 




12(12) 






4(2) 
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TABLE IV. Distribution of topological charge along with the prediction based on a 
non-interacting gas of instantons and anti-instantons using the average from the Poisson distri- 
bution. 





SU(2) 


SU(2) 


SU(2) 


SU(2) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


SU(3) 


a 


2.3 


2.4 


2.4 


2.5 


5.71 


5.75 


5.75 


5.75 


5.85 


5.85 


L 


8 


8 


16 


16 


16 


8 


12 


16 


8 


16 


N 


50 


200 


200 


200 


200 


400 


200 


400 


400 


206 


Q = 


19(17) 


75(76) 


61(62) 


134(136) 


38(36) 


293(297) 


103(96) 


101(110) 


379(378) 


153(143) 


Q = l 


21(22) 


23(22) 


82(86) 


60(57) 


67(65) 


101(95) 


72(82) 


164(166) 


19(22) 


44(57) 


Q = 2 


7(8) 


2(2) 


41(37) 


6(6) 


45(47) 


6(8) 


22(19) 


86(83) 


2(0) 


9(6) 


g = 3 


2(2) 




14(11) 




26(28) 




3(3) 


35(30) 






Q = 4 


1(0) 




2(3) 




18(14) 






12(9) 






g = 5 










3(6) 






2(2) 






g = 6 










1(2) 












g = 7 










1(1) 












g = 8 






















g = 9 






















Q = 10 










1(0) 












{Q')/v 


1.66 


0.31 


0.25 


0.05 


0.64 


0.31 


0.28 


0.33 


0.07 


0.05 



Finally, the assumption of a non-interacting gas of instantons and anti-instantons leads 
to a simple expression for the distribution of topological charge, 



T(Q)=e-<")/Q((n)) 



(14) 



where Iq is the modified Bessel function of order Q. Our data can be cross-checked against 
this prediction also, and this is shown in Table 0. Again, we see fairly good agreement 
with the data. We also note that (Q^) = [n) under this assumption, and this is essentially 
the case when we compare Table y and Table |rV|. 

This simple model of associating the zero and the small, non-zero eigenvalues with inde- 
pendent instantons and anti-instantons seems to account for the main statistical properties 
of the n+,n_ data. 



IV. SMALL EIGENVALUES AND A CHIRAL CONDENSATE 

This section deals with the spectrum of the small, non-zero eigenvalues. We begin with 
a discussion relating the spectrum to expectations from the instanton — anti-instanton gas 
model. This is followed by an analysis of the way in which these eigenvalues contribute to 
physical quantities such as the chiral condensate. 

In the previous section, we were only concerned with the number of small eigenvalues 
(eigenvalues below 0.05) in each configuration. Our computation of the low lying eigenval- 
ues of H^ along with the spectral flow of H^^ enables us to obtain this number with very 
good certainty per configuration. Because of the numerical accuracies involved, the small 
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eigenvalues can be trusted only with an absolute accuracy of 0.005, and we can create a 
histogram of the low lying eigenvalues with a bin width of 0.005 giving us a total of ten bins 
in each ensemble. 

A detailed analysis of the relationship of small eigenvalues, instantons, and the conden- 
sate at zero temperature can be found in [^. It was shown there that the spectrum from a 



model with fermions in a background of dilute instantons and anti-instantons gives a chiral 
condensate and is consistent with chiral random matrix theory. Some discussion of the finite 
temperature spectrum appears in Pj. We have used a much simpler version of these ideas. 
The most important effect for widely separated instantons and anti-instantons is from 
the mixing of the n would-be zero modes of the fermions. At infinite separation, there would 
be n modes with A = 0. At finite temperature and in the continuum, the spatial tails of these 
wave functions are exponential. Further a configuration with e.g. n = n+ has that number 
of exact zero modes. This leads us to a model with instantons and anti-instantons that are 
still distributed independently over the configurations and with an n x n interaction matrix 
for the n modes that is made up of n_|_ x n+ and n_ x n_ blocks on the diagonal and n_ x n+ 
and 72+ X n^ off-diagonal blocks. The diagonal blocks corresponding to instanton — instanton 
fermion mode interactions and to anti-instanton — anti-instanton fermion mode interactions, 
including self-interactions are zero. 

Tij = if z and j are both instantons or both anti-instantons. (15) 

In the "off-diagonal" blocks there is an entry for each ordered way that an instanton can be 
paired with an anti-instanton. It has the form 

Tij = hoe^'^^'^'^''^ if {i,j) is an instanton — anti-instanton pair. (16) 

The energy scale is determined by the constant ho. The distance between the instanton — 
anti-instanton pair is d{i,j). The length scale of the mode interactions is D. We expect D 
to be of order A^^ since that is the range of the continuum mode tails. 

The next step is to generate data from the model and compare it with the real data. The 
model calculation is in the continuum with periodic boundary conditions, so that d{i,j) is 
the shortest distance connecting the pair. The {n^,n_) values are distributed according to 
Eqn. ( p!2D and Eqn. ([T3| ) with (n) obtained from Table |1[ The positions of the instantons 
and anti-instantons are selected at random. Then the matrix elements of T are computed. 
Finally the eigenvalues of T are computed. Note that for a matrix with the form of T, 
the spectrum has \Q\ zeros and (n — \Q\)/2 pairs ±A. This process is repeated for many 
configurations. The resulting set of eigenvalue distributions is then compared with the data 
from the various ensembles. The conclusion is that a value oi D = 2 best fits the data 
for all the ensembles. We have shown the comparison between the "real" data and the toy 
model for two of the cases in Fig. 0. Thus the spectrum of the small eigenvalues can also be 
explained by the simple instanton — anti-instanton gas picture. 

A physical role for these modes is suggested by the fact that the chiral condensate is 
determined by the infinite- volume fermion spectral density at zero eigenvalue. The most 
straightforward thing to do would be to use the computed eigenvalues to do the sums in 
Eqns. (1^) - (|Tl]) for the condensate, a;, and Xao- However, we need the /x — *> limit, and it 
is apparent from the forms of the sums that they are very sensitive to the positions of the 
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FIG. 2. Histogram of the small eigenvalues along with the histogram from the toy model for 
£' = 1,2,3 with the appropriate value of (n) from Table || for the SU(2) ensemble at /? = 2.4 and 
SU(3) ensemble at /3 = 5.75. Both lattices have a spatial volume of 16'^. 
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eigenvalues when /x is very small. As mentioned earlier in this section, the small eigenvalues 
are not precisely known, so that the value of the sum is unreliable for fi less than about 
0.01, and this is not small enough to be useful. 

However, it is possible to compute the chiral condensate, u, and Xao by a stochastic 
method using the overlap Dirac operator. The method is specific to the overlap Dirac oper- 
ator and is described in Ref . [|15,0 . Using this method, we have computed the condensate, 
the anomalous U(1)a breaking, and the scalar susceptibility for the SU(2) ensemble with 
(3 = 2.4 on the 16^ lattice and for the SU(3) ensemble with j3 = 5.75 also on the 16^ lattice. 
The results are shown in Fig. ^. If the steep decrease of the condensate toward zero below 
H = 0.002 is interpreted as a finite volume effect, then these data are evidence for chiral 
symmetry breaking. Additional evidence for chiral symmetry breaking can be seen in the 
plots for u. We see there is a quark mass region in the non-topological contribution to 
uj from /i = 0.002 to 0.05 where the expected divergence of - sets in. The region below 
/i = 0.002 is where the finite volume effects become large. 

The shapes of the curves in Fig. |^ are consistent with the qualitative features of the 
computed eigenvalues: small eigenvalues concentrated close to A = 0, a sparsely populated 
region around A = 0.05, and the dense eigenvalues of the bulk of the spectrum beginning 
at larger eigenvalue. The strongest effect is a dip in Xao due to the presence of two terms 
in the sum in Eqn. (^) with opposite signs along with a spectral distribution that has two 
regions (a small eigenvalues region and a bulk region) separated by a sparsely populated 
region around A = 0.05 as seen in Fig. [I|. 

One may ask whether the simple model, which reproduces the properties of the small 
eigenvalues, gives a chiral condensate. The short answer to that question is "Yes". Our 
model is just a simplification of the model in Ref. pSf , which does give a condensate. The 
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longer answer comes after testing the simplified model itself for the necessary properties. 
With {n)/V fixed and increasing V, the staircase function 

N{\V)= f\ap{a,V) (17) 

in the large V and small A region needs to approach a function of AV^ that is linear for large 
W . N{X, V) as computed with modest statistics from the model eigenvalues is consistent 
with that behavior. Thus it appears that the eigenvalues in the "small region" of the 
spectrum are sufficient to give rise to a condensate. 



V. CONCLUSIONS 

We have studied the spectrum of the overlap Dirac operator on several quenched ensem- 
bles in the deconfined phase. We find clear evidence for topology in the deconfined phase. 
We also found that the topological susceptibility decreases sharply as we go deeper into the 
deconfined phase. The gauge field contains instanton and anti-instanton like objects which 
are well described as a dilute gas. This is supported by a part of the spectrum of the overlap 
Dirac operator consisting of "small" (in our case < 0.05) eigenvalues well separated from 
the bulk. This small part is consistently described as arising from a dilute gas of instantons 
and anti-instantons. It appears that the small eigenvalues produced by such a dilute gas 
are sufficient to create a chiral condensate and hence spontaneous chiral symmetry breaking 
even in the deconfined phase of quenched gauge theories. 
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FIG. 3. The stochastic estimates of {ipip), uj and Xa^ for the SU(2) ensemble at (3 = 2.4 and 
SU(3) ensemble at /3 = 5.75. Both lattices have a spatial volume of 16^. 
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